Compression of ultrashort UV pulses in a self-defocusing gas 
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Compression of UV femtosecond laser pulses focused into a gas cell filled with xenon is reported 
numerically. With a large negative Kerr index and normal dispersion, xenon promotes temporal 
modulational instability (MI) which can be monitored to shorten ~ 100 fs pulses to robust, singly- 
peaked waveforms exhibiting a fourfold compression factor. Combining standard MI theory with a 
variational approach allows us to predict the beam parameters suitable for efficient compression. At 
powers < 30 MW, nonlinear dispersion is shown to shift the pulse temporal profile to the rear zone. 

PACS numbers: 42.65.Tg, 42.65.-k,52.38.Hb,42.68.Ay 



Since the nineties, considerable progress have been re- 
ported in the field of optical pulse compression [l[ . From 
the first achievements of few-cycle pulses [H, the om- 
nipresent idea has been to control spectral broadening, 
in order to reach shorter light structures. This goal can 
be completed by post-compressing pulses resulting from 
the interplay between self-phase modulation (SPM) and 
group-velocity dispersion (GVD) in hollow waveguides 
and soliton compression devices, e.g., fibers. Other tech- 
niques have been proposed, using negative phase shifts 
based on cascaded quadratic nonlinearities, which can 
preserve most of the pulse energy [sl, 0, [5[. Alterna- 
tively, few-cycle filaments of light can be created from 
the balance between Kerr self-focusing and plasma defo- 
cusing at higher peak powers 0, 0, S] • Frequency conver- 
sion processes, that enlarge more the spectrum, have also 
been exploited to generate ultrashort pulses at UV wave- 
lengths. Recently, Fuji et al. 0] produced 12-fs pulses 
at 260 nm with ~ 10 /iJ energy through the four- wave 
mixing of 400 and 800 nm pulses in filamentation regime 
13, El. The same scheme was numerically optimized to 
compress pulses below 2 fs in low-pressure argon cell [HI . 

Apart from the effective self-defocusing nonlinearity 
promoted by cascaded quadratic interactions few at- 
tention has been paid to materials having a negative Kerr 
response. These exist however, such as helium and xenon, 
for laser wavelengths in the range of ~ 238 — 249 nm. In 
[l3| , KrF laser light was used to measure the Kerr index 
712 of xenon, which can attain important negative values 
close to a two-photon resonance at 249.6 nm. In [3], nu- 
merical evidence was given to the principle of pulse com- 
pression in (1+1) dimensional hollow fibers filled with 
xenon, where pulses at 243 nm could be drastically short- 
ened with minimum energy losses. In higher dimensional 
systems, nonlinear defocusing can offer a rich variety of 
dynamical patterns, e.g., ring dark solitons and vortices 
[l5[. For a medium supporting normal GVD, modula- 
tional instability (MI) can moreover develop in time, pro- 
viding further potentiality for pulse compression. 

In this Letter, we investigate new propagation regimes 



implying a defocusing nonlinearity without any guiding 
device. 100 fs Gaussian pulses at 243 nm central wave- 
length can be compressed in a cell of xenon at ambi- 
ent pressure to about 25 fs along a stage of modula- 
tional instability. These compressed structures support 
propagation ranges up to 1 m. Because nonlinearity is 
highly dispersive in the UV, we take into account the 
full frequency dependency of the third-order susceptibil- 
ity x^^\-^]^"' where uo = uo'" uo" uo' . Re- 
sults obtained from such a rigorous approach are still 
compatible with those inferred from a classical nonlinear 
Schrodinger (NLS) model. Qualitative behaviors can be 
predicted by combining a variational method with 
the classical MI theory for plane waves l3, [13] • 

We assume a linearly-polarized electric field ~ 
g^zkoz-tujot _^ where E{r^z^t) denotes its envelope 
and / = is the pulse intensity. The wave number 
^0 = nocoo/c involves the linear refractive index of the 
gas, no — 1, at central frequency coq. Because the max- 
imum intensity /max never exceeds 0.5 TW/cm^ in the 
coming simulations, the peak density of free electrons 
created by ionization always remains below 10^^ cm~^, 
rendering plasma generation negligible. The physics can 
thus be described by the following envelope equation in 
Fourier space 
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where r = \f x- 



^ + vp" and a; = — cjo denotes the enve- 
lope frequency. Here, z is the propagation variable and t 
a retarded time. Linear dispersion for xenon is included 
via k{uj) [13], while the expression for the nonlinear sus- 
ceptibility x*^^^ is that given in [IJ] [see also Fig. [41(b)]. 
For sufficiently narrow spectral bandwidths, nonlinear as 
well as high-order linear dispersion, self-steepening and 
space-time focusing operators are expected to have a mi- 
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nor influence on the qualitative pulse dynamics. We will 
thus compare results obtained from Eq. ([T]) with those 
from the classical NLS equation 
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with = 13.12 fsVcm [19] and = -1.56 x 

10~^^ cm^/W Both propagation models are inte- 

grated numerically for Gaussian pulses with input power 
Pin, beam waist wq and 1/e^ pulse half- width Tp. The 
delocalizing action of the nonlinearity is balanced by lin- 
early focusing the incident beam. For notational con- 
venience, we will employ the deflnition of the critical 
power for the self- focusing of Gaussian beams, Per — 
AQ/27rno|n2| — 6 MW. Key processes should thus be 
transverse diffraction, normal GVD and Kerr defocus- 
ing. This allows us to capture the nonlinear dynamics by 
means of two different, simple analytical techniques. 

On the one hand, we perform a two-scale variational 
analysis fid\ resulting into the dynamical system 
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where p = Pin /Per- This system governs the beam waist 
[w{z)] and 1/e^ temporal radius [r{z)] of Gaussian pulses, 
starting from w{0) = wq and r(0) = Tp. The normalized 
on-axis intensity behaves as I / Iq = WQrp/[w'^{z)T{z)] and 
the lens action can be modelled through suitable phase 
contributions containing dzw{0) = —wo/f. Such an 
approximation method cannot describe the flne spatio- 
temporal deformations of the pulse. Nevertheless, it usu- 
ally yields estimates of the maximum intensity and of 
the pulse scales in space and time, which support the 
comparison with direct numerical results. As seen from 
Eq. (|3]), the pulse dynamics differs from standard self- 
focusing by the sign in front of the nonlinear term ~ p. 
When the pulse duration does not vary too much, e. g., 
for weak nonlinear it ies, the equation for w{z) can readily 
be integrated and yields a focus reached at the distance 
^min - //[I + (//2;o)^(l +p/V^)], zo = TTUow^ / Xq being 
the Rayleigh length of the input beam. At high enough 
powers, the pulse can, in contrast, undergo significant 
compression in time, as the Kerr response competes with 
normal GVD. This property will be exploited below. 

On the other hand, the standard stability analysis for 
plane waves can bring further insight into the pulse dy- 
namics. We assume that the highest intensity zones 
of the pulse serve as plane- wave distributions, as long 
as local perturbations have typical wavelengths (peri- 
ods) less than the size (duration) of the background 
field with intensity close to its maximum, /max- We 
can thus linearize Eq. ([2]) for perturbations oscillat- 
ing with the transverse wavenumber k± and frequency 
uj [13, so that the MI growth rate expresses as 7 = 
Re (n ^2^0 1 n2 \ /max /c -n^), where Q'^ = k^'^^uj'^ /2 - 



k\/2kQ. Modulational instability develops for positive 
values of only. Following [13], a necessary condition 
for MI is that the optimum perturbation wavenumber 
/c^^^ and frequency aJmax, linked to each other by 
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must satisfy k^^^ > \/2i\lw^ and cJmax > V2ti /Tp at 
given /max- MI fully takes place whenever perturba- 
tion frequencies uJ > cJmax exist. Fixing the value of 
^max ^ \/27r/i(;min , whcrc Wmin is the smallest beam 
waist reached at focus, it is then sufficient to evaluate 
whether the frequency range \/27r/rp < uJ < cJup = 
\/27r/rmin fulfllls the condition uJ^p > cJmax, to conclude 
if MI is efficiently seeded. Here, cJup includes the min- 
imum pulse duration Tmin- We conjecture that reliable 
estimates for /max, '^min and Tmin are provided by the 
variational equation (|3]). In the following, these quanti- 
ties will be used to evaluate /c^^^, cJup and cJmax- 

An example for a stable conflguration is shown in 
Fig. [U where a Gaussian pulse with Pin = 6 GW, 
Wo = 1 mm and Tp = 200 fs is simulated. From our varia- 
tional model we find /max — 40 GW/cm^, Wmin — 100 /im 
and Tmin — 100 fs, resulting in k^^^ 444 cm~^, 
oJup — 0.044 fs~^ and oJmax — 0.29 fs~^. Since oJup < 
^max, we expect stability of the temporal stripe designed 
by the pulse in the (t, z) plane. The peak intensity 
[Fig. HJa)], computed from direct simulations (solid and 
dashed curves) and from Eq. ([3]) (dotted curve), reaches a 
maximum near the focal distance 2: = 1 m. Here, the in- 
tensity curves accounting or not for nonlinear dispersion 
are almost identical. Figure [T](b) shows the spatial profile 
computed from Eq. ([1]). The beam focuses then diffracts 
with a divergence angle given by tan^beam ^ ^o/7TWmin- 
Figure [TJc) illustrates the evolution of the on-axis tempo- 
ral profile in the (t, z) plane, which confirms robustness 
of the temporal profile against MI. The white dots re- 
produce the functions w{z)^ t{z) computed from Eq. ([3]). 
In Fig. [TJd), the minimal pulse duration of about 100 fs 
is reached after a propagation distance of 4 m for both 
Eqs. ([1]) and ([2]), while the same extent in time is attained 
from the variational model at z = 8 m. 

Let us then inspect regions of temporal instability. 
When we increase both the input power and focal length 
by a factor five, the variational model predicts /max — 
6 GW/cm^ 

, '^min — 640 /im and Tmin — 33 fs, resulting in 

^max _ ^^-1^ ^ q fg-1 and cJmax ^ 0.071 fs-^ 

Since cJup > ^max, nonlinearities should thus seed MI and 
compression in time. This is actually confirmed by the 
simulations shown in Fig. [21 Compression leads to a short 
peak of ~ 25 fs duration at z = 4 m, before MI fully 
breaks up the pulse temporal distribution. The modu- 
lation period is about 27r/cJmax ~ 150 fs, which agrees 
with Figs. Efd) and (e). Note that again the minimum 
duration predicted by Eq. (j3|) occurs beyond the com- 
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FIG. 1: (a) Peak intensities of a Gaussian pulse with Pin = 
6 MW, Wo = 1 mm and Tp = 200 fs, focused in xenon at atmo- 
spheric pressure (/ = 1 m). Sohd hne corresponds to the com- 
plete model Eq. ([T]), dashed (superimposed) line to Eq. (|2]), 
and dotted line to the variational model Eq. (|3)). (b) Spatial 
and (c) temporal dynamics of the full model. White dots re- 
produce the beam waist and 1/e^ temporal extent computed 
from Eq. (|3|). (d) Power profiles after 4 m of propagation 
for the full model (solid line) and the reduced model Eq. (|2]) 
(dashed line). The blue line shows the initial pulse. 

pression stage seen in the simulations. The numerical 
pulse is subject to leaks of power, which the variational 
approach cannot describe by preserving the input power 
within a Gaussian ansatz [l6|. However, Eq. (|3]) repro- 
duces compression rates agreeing well with the numer- 
ical data. While nonlinear dispersion clearly influences 
the pulse dynamics by moving its temporal centroid to- 
wards positive times, we observe qualitative agreement 
between results obtained from Eqs. ([1]) and ([2]) [compare 
Figs.Hb), (c) andUd), (e)]. 

Figure [3] concerns the main result of this Letter. By 
tuning the input parameters appropriately, MI can be 
handled, such that only one modulation mainly affects 
the pulse temporal profile. Since a negative n2 may lead 
to compression at distances z > Zmin, the pulse can then 
not only reach very short durations, but also keep them 
over long distances. The condition for this regime of op- 
timal pulse compression is obviously cJup — cJmax- Let 
us examine the evolution of a Gaussian pulse, for which 
we reduce both initial waist, pulse duration and power 
to Wo = 0.3 mm, Tp = 100 fs, and Pin = 24 MW, respec- 
tively, to obtain cJmax — 0.14 fs~^ and oJup — 0.17 fs~^. 
Figures [D^a) and (b) detail the pulse evolution computed 
from Eqs. ([2]) and ([T]), respectively. MI is weakly seeded 
and the 100 fs pulse undergoes three modulations, among 
which only the central one is amplified and compressed 
along further propagation down to ~ 25 fs. The resulting 
structure appears to behave like a ID sech-soliton in time, 
capable of preserving high intensity and short durations 
over at least 50 cm. Again we observe a slight shift to 
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FIG. 2: (a) Peak intensities of the same Gaussian pulse as in 
Fig. [T] but with Pin = 30 MW and / = 5 m, using the same 
plotstyle as in Fig. [T]^a). (b) Temporal on-axis dynamics ob- 
tained from Eq. (|2|; white dots correspond to the variational 
approximation, (c) Same from the complete model Eq. ([T]). 
Spatio-temporal intensity distribution at 2; = 8 m showing MI 
from (d) the reduced model, (e) the complete model. 



the trailing part of the pulse caused by the nonlinear dis- 
persion. At such weak intensities, we can also infer that 
two-photon absorption remains of small influence. It is 
nevertheless important to point out that, in contrast to 
filamentary compression in self-focusing regime [6, 7,0], 
short pulses obtained with the above compression scheme 
are homogeneous in radial direction. Hence, compression 
applies not only to the on-axis intensity profiles, but also 
to the pulse power profiles. Moreover, the XFROG trace 
shown in Fig. [SJd) indicates that the compressed pulses 
are nearly transform limited. 

A last point concerns spectral broadening, as 
SPM is responsible for frequency variations Aco :^ 
/coAz|n2|/max/Tmin arouud the focus. Spectral variations 
should remain confined within the narrow bandwidth of 
238-249 nm (7.57-7.92 PHz), in order to preserve a neg- 
ative value of n2 ~ x^^H^5 ^5 ^) along the optical 
path [see Fig. [4](b)]. Figure Ht^a) confirms that the spec- 
trum has a small bandwidth of Aco < 0.1 PHz at half- 
maximum. For the simulation shown in Fig. [3l spectral 
wings reach the range where 112 is not necessarily neg- 
ative, but the dominant peak always remains located 
around the central frequency [solid line in Fig. (H^a)]. 
We can notice that the usual broad low-intensity super- 
continuum (Accuse /^o > 1) routinely observed from the 
propagation of plasma-induced filaments and caused by 
the mechanism of intensity clamping is missing. We in- 
deed find AcJsc/^o < 0.1, where Acjsc is measured at 
10~^ times the maximum spectral intensity. Supercon- 
tinuum generation is thus inhibited, as the peak intensity 
stays at moderate values. 

Further increase of the input power to 60 MW leads 
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simple analytical procedures, we proved that this new 
compression mechanism follows from modulational in- 
stability in time, which can be controlled to optimize 
the self-compression process over long distances. To end 
with, we find it worth emphasizing that nonlinear disper- 
sion yielding variations in the Kerr index of at least one 
order of magnitude does not significantly change results 
compared to propagation models assuming a constant 
susceptibility. Here, we indeed showed that at moderate 
powers deviations associated with nonlinear dispersion 
remain limited and their qualitative effect is similar to 
classical pulse self-steepening in focusing media. 

This work has been performed using HPC resources 
from GENCI-CCRT/CINES (Grant 2009-x2009106003). 



FIG. 3: Dynamics of a Gaussian pulse with Pin = 24 MW, 
wq = 0.3 mm and Tp — 100 fs, focused in xenon with / = 1 m. 
(a) Temporal on- axis dynamics obtained from the reduced 
model; white dots correspond to variational results, (b) Same 
from the full model, (c) Power profiles after 1.5 m of propa- 
gation for the full model (solid line) and the reduced model 
(dashed line). The blue line shows the initial pulse, (d) 
XFROG trace at 2; = 1.5 m using a 25 fs reference pulse. 
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FIG. 4: (a) On-axis normalized spectra for the pulses shown 
in Fig.^e) (dashed line) and Fig.[3i;c) (solid line). The dotted 
line represents the spectrum of the 100 fs pulse ai z = 0. (b) 
Real part of the normalized nonlinear susceptibility X<3> in 
the degenerated case (see Ref. [3). 



to the occurrence of strong MI and multiple pulse split- 
ting (cJup > ^max, not showu here). For such high pow- 
ers, however, a major part of the broadened spectrum 
falls into the range of multi-photon resonances below 
7.57 PHz, and strong discrepancies between the solutions 
of Eq. ([1]) and Eq. ([2]) appear. Consequently, while it is 
possible in the NLS model to produce even shorter pulses 
< 20 fs by simply increasing the initial peak power, dis- 
persion of the nonlinear susceptibility prevented us from 
achieving shorter durations with Eq. ([1]). Such propaga- 
tion regimes will be addressed in a future publication. 

In summary, numerical simulations have highlighted 
the dynamics of UV pulses focused into a self-defocusing 
gas cell. Ultrashort optical structures can naturally be 
formed in (3+l)-dimensional media with negative 712 and 
normal dispersion, and propagate over long ranges with- 
out experiencing a wide spectral broadening. By mixing 
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